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ABSTRACT 

Useful  constraints  on  the  orbits  and  mass  ratios  of  astrometric  binaries  in  the  Hipparcos  catalog  are  derived  from 
the  measured  proper  motion  differences  of  Hipparcos  and  Tycho-2  (A/k),  accelerations  of  proper  motions  ( jj),  and 
second  derivatives  of  proper  motions  (/2).  It  is  shown  how,  in  some  cases,  statistical  bounds  can  be  estimated  for  the 
masses  of  the  secondary  components.  Two  catalogs  of  astrometric  binaries  are  generated,  one  of  binaries  with 
significant  proper  motion  differences  and  the  other  of  binaries  with  significant  accelerations  of  their  proper  motions. 
Mathematical  relations  between  the  astrometric  observables  A/i,  fi,  and  /i  and  the  orbital  elements  are  derived  in  the 
appendices.  We  find  a  remarkable  difference  between  the  distribution  of  spectral  types  of  stars  with  large  accel¬ 
erations  but  small  proper  motion  differences  and  that  of  stars  with  large  proper  motion  differences  but  insignificant 
accelerations.  The  spectral  type  distribution  for  the  former  sample  of  binaries  is  the  same  as  the  general  distribution 
of  all  stars  in  the  Hipparcos  catalog,  whereas  the  latter  sample  is  clearly  dominated  by  solar-type  stars,  with  an 
obvious  dearth  of  blue  stars.  We  point  out  that  the  latter  set  includes  mostly  binaries  with  long  periods  (longer  than 
about  6  yr). 
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1.  INTRODUCTION 

Binaries  with  invisible  or  unresolved  companions  may  create 
problems  in  the  data  reduction  of  a  large  and  accurate  astro- 
metric  catalog.  Owing  to  the  orbital  motion,  which  may  be 
measurable  astrometrically,  these  stars  cannot  be  accurately 
described  by  the  five-parameter  model  that  includes  two  posi¬ 
tion  components,  two  components  of  proper  motion,  and  par¬ 
allax  as  unknowns.  Instead,  a  more  complex  model  of  seven  or 
nine  free  parameters  (including  the  acceleration  fi  and  second 
derivative  of  proper  motion  fl)  or  a  complete  orbital  adjustment 
of  12  parameters  may  be  necessary.  About  2.2%  of  stars  in  the 
Hipparcos  catalog  require  such  a  special  treatment  (Perryman 
et  al.  1997).  With  the  anticipated  greater  accuracy  of  future  as¬ 
trometric  missions  such  as  the  Space  Interferometry  Mission 
(SIM)  and  Gaia,  the  fraction  of  stars  exhibiting  nonlinear  ap¬ 
parent  motions  will  dramatically  increase  (Kaplan  &  Makarov 
2003). 

While  astrometric  ramifications  of  the  problem  have  been 
addressed  in  a  number  of  publications,  less  attention  has  been 
paid  to  the  prospect  of  characterizing  new  binaries  through  these 
astrometric  effects.  A  catalog  of  about  7000  A/r  binaries — those 
with  significant  differences  in  the  Hipparcos  and  Fifth  Funda¬ 
mental  Catalogue  proper  motions — was  created  by  Wielen  et  al. 
(2000).  Gontcharov  et  al.  (2001)  produced  a  similar  catalog  of 
a  few  hundred  astrometric  binaries  with  nonlinear  motion.  In 
further  development  of  the  approach  used  in  these  previous  in¬ 
vestigations,  we  derive  formulae  to  estimate  the  physical  size  of 
the  orbit  (i.e.,  the  semimajor  axis  in  AU)  for  astrometric  binaries 
with  invisible  companions  from  the  available  astrometric  data  in 
the  Hipparcos  and  Tycho-2  (Hog  et  al.  2000)  catalogs.  The  bi¬ 
naries  are  divided  into  three  categories:  (1)  stars  with  statistically 
significant  differences  between  the  short-term  Hipparcos  proper 
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motions  and  the  long-term  Tycho-2  proper  motions;  (2)  stars 
with  measured  acceleration  of  proper  motions  as  given  in  the 
Hipparcos  catalog.  Annex  DMSA/G;  and  (3)  stars  with  mea¬ 
sured  second  derivatives  of  proper  motion  (besides  the  first  de¬ 
rivative,  or  acceleration)  as  given  in  the  same  annex.  We  present 
two  catalogs  of  astrometric  binaries,  one  of  which  includes  A p 
binaries,  and  the  other  fi  and  jl  binaries. 

Our  catalogs,  which  are  mostly  compilations  of  previous  as¬ 
trometric  data,  are  intended  for  researchers  interested  in  binary 
stars,  in  particular  for  the  search  of  binary  brown  dwarfs  and 
other  low-mass  companions.  The  catalogs  include  only  reliably 
identified  long-period  astrometric  binaries,  which  are  difficult  to 
detect  by  other  means.  The  objectives  of  the  previous  Heidelberg 
catalogs  TYC2+HIP  and  ARIH1P  (Wielen  et  al.  2001a,  2001b), 
on  the  contrary,  were  ( 1 )  to  provide  the  most  accurate  proper  mo¬ 
tions  for  certain  samples  of  stars  and  (2)  to  obtain  “mean”  proper 
motions.  The  emphasis  of  the  ARIH1P  catalogs  was  heavily  on 
the  “astrometrically  excellent”  stars  that  could  be  used  as  an 
optical  representation  of  the  International  Celestial  Reference 
System.  Our  criteria  are  strict  in  the  selection  of  high-fidelity 
binary  stars,  whereas  the  ARIHIP  criteria  select  nonbinary 
“excellent”  stars.  (Essentially,  for  identifying  binary  systems, 
the  difference  is  between  minimizing  false  positives  and  mini¬ 
mizing  false  negatives.)  Most  importantly,  the  Heidelberg  team 
empirically  estimated  an  all-sky  correlation  between  the  Hip¬ 
parcos  catalog  and  Tycho-2  catalog  proper  motions.  Since  this  is 
a  positive  correlation,  applying  it  to  the  detection  of  stars  with 
perfectly  unperturbed  motion  makes  it  easier  to  achieve  the  sta¬ 
tistical  binarity  threshold.  The  real  correlation  may  differ  from 
the  mean  value  at  different  parts  of  the  sky  because  of  varying 
number  density,  brightness  of  the  Hipparcos  stars,  and  system¬ 
atic  errors  (Wielen  et  al.  2001b).  While  this  approach  may  be 
justified  if  one  wants  a  reliable  sample  of  single  stars,  the  filtered 
perturbed  stars  are  not  necessarily  binaries,  and  certainly  not 
necessarily  the  astrometric  binaries  with  faint  companions  that 
are  the  focus  of  this  paper. 

Parameters  Qq  (§  2)  and  Q\  (§  3)  are  computed  for  all  stars  in 
our  catalogs.  They  are  related  to  the  lower  bounds  of  the  mass 
ratio,  and  their  use  is  exemplified  in  §  2  by  the  young  binary 
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TABLE  1 

A  Catalog  of  A/t  Binaries 


Hipparcos  Tycho 

-  -  A 

HIP  No.  TYC1  TYC2  TYC3  ua »  ns  ua,  n6  n  Binary?  logQ0 

(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10)  (11) 


68 .  1178  1142  1  -99.7  ±1.1  -315.9  ±0.8  -69.0  ±  0.9  -305.0  ±  0.9  31.8  ±1.2  A  -0.8 

93 .  4663  456  1  54.9  ±1.3  -73.3  ±  0.7  50.8  ±  0.9  -68.7  ±  0.9  16.7  ±1.3  -1.2 

171 .  1732  2731  1  778.6  ±  2.8  -918.7  ±  1.8  829.9  ±  1.2  -989.4  ±1.1  80.6  ±  3.0  -0.8 

290 .  6989  1159  1  93.7  ±1.1  -42.8  ±  0.6  92.5  ±1.4  -56.7  ±1.4  15.3  ±  1.0  -0.8 

305 .  6418  1218  1  88.6  ±1.1  -17.7  ±0.9  100.4  ±1.8  -20.3  ±  1.4  20.4  ±1.0  A  -1.0 

329 .  5263  680  1  143.7  ±  1.9  66.6  ±  0.8  149.9  ±  1.7  75.2  ±1.4  9.2  ±1.5  -0.7 

356 .  4669  211  1  -42.5  ±1.2  -17.2  ±0.6  -33.8  ±  1.1  -20.3  ±1.2  9.9  ±1.1  -0.8 

359 .  9346  1115  1  177.5  ±  0.9  -15.2  ±0.7  161.6  ±1.5  -16.5  ±1.4  17.1  ±0.9  -0.8 

457 .  7529  512  1  -12.8  ±0.8  —1.3  ±  0.7  -14.0  ±0.9  3.0  ±  0.9  2.0  ±1.0  -0.5 


Notes. — Several  entries  are  shown  as  an  example;  the  complete  version  of  the  catalog  is  available  through  the  CDS  Web  site.  Col.  (1):  Hipparcos  numbers.  Cols.  (2)- 
(4):  Tycho  identification  numbers  TYC1,  TYC2,  and  TYC3.  Cols.  (5)  and  (6):  Hipparcos  proper  motion  components  and  p,f  in  mas  yr-1  and  their  standard  errors. 
Cols.  (7)  and  (8):  Tycho-2  proper  motion  components  and  their  errors.  Col.  (9):  Hipparcos  parallaxes  and  their  errors  in  mas.  Col.  ( 10):  Flag  “A”  is  set  for  stars  that  are 
also  //  binaries.  Col.  (11):  log  Qo  factors  for  the  lower  limit  of  q2- 


AB  Dor  with  a  brown  dwarf  companion.  On  the  contrary,  the 
“cosmic  error”  given  in  the  Heidelberg  compilations  is  a  measure 
of  astrometric  disturbance  and  is  not  related  directly  to  the  orbital 
elements.  In  the  appendices,  we  derive  formulae  linking  the  ob¬ 
servables  Ap,  p,  and  p  with  the  orbital  parameters.  This  anal¬ 
ysis  will  be  useful  in  astrometric  applications.  For  example,  a 
number  of  SIM  target  stars  and  grid  stars  will  have  to  be  reduced 
with  the  seven-parameter  astrometric  model  that  includes  accel¬ 
eration  components,  and  the  formalism  is  directly  applicable  in 
the  ensuing  astrometric  analysis. 

2.  Ap  BINARIES 

Binarity  of  this  type  is  revealed  when  two  astrometric  catalogs 
of  proper  motion  are  compared.  Depending  on  the  accuracy  of 
the  catalogs,  a  certain  fraction  of  stars  will  show  differences  in 
proper  motion  in  excess  of  a  statistically  plausible  measurement 
error.  This  happens  when  one  of  the  catalogs  includes  short-term 
proper  motions  (i.e.,  based  on  observations  collected  during  a 
relatively  short  time),  while  the  other  is  based  on  long-term  ob¬ 
servations  of  star  positions.  In  binaries  of  sufficiently  long  pe¬ 
riods,  the  reflex  orbital  motion  of  the  primary  will  be  captured  in 
the  observed  short-term  proper  motions.  The  long-term  proper 
motion,  on  the  other  hand,  will  be  closer  to  the  true  center-of- 
mass  motion  of  the  system.  Thus,  we  define  A p  binaries  as  stars 
that  have  instantaneous  (or  short-term)  proper  motions  signifi¬ 
cantly  different  from  the  quasi-inertial  motion  of  the  center  of 
mass.  For  the  majority  of  stars,  the  unknown  astrometric  excur¬ 
sion  of  the  photocenter  from  the  center-of-mass  position  will 
be  averaged  out  to  quantities  comparable  to  or  smaller  than  the 
astrometric  formal  error.  Because  of  the  averaging  effect  in  long¬ 
term  proper  motions,  it  is  expected  that  mostly  long-period  bi¬ 
naries  are  detectable  as  Ap  binaries.  Some  short-period  binaries 
can  accidentally  produce  significant  proper  motion  difference 
too  if  their  apparent  astrometric  excursions  are  large  and  only  a 
few  observations  are  available. 

Our  catalog  of  A p  binaries  ( Table  1 )  includes  1 929  stars  with 
Ap  differences  between  the  Hipparcos  and  Tycho-2  catalogs 
larger  than  3.5  a  in  at  least  one  of  the  coordinate  components. 
A  typical  observational  time  span  for  a  Hipparcos  star  is  about 
3.2  yr,  while  Tycho-2  proper  motions  are  derived  from  astro- 
metric  observations  often  spanning  over  a  century  (Urban  et  al. 
2000).  Since  Hipparcos  proper  motions  were  used  for  correct¬ 
ing  the  system  of  ground-based  catalogs,  they  are  positively 


correlated  with  the  Tycho-2  proper  motions  for  common  stars. 
These  correlations  are  not  known  for  each  star,  and  we  chose 
to  neglect  them  in  order  to  reduce  the  probability  of  false  pos¬ 
itives,  that  is,  to  generate  a  more  reliable  set  of  binary  stars. 
Because  of  the  neglected  positive  correlation,  the  calculated 
Ap  is  smaller  than  the  true  value,  which  tends  to  give  an  un¬ 
derestimate  of  the  secondary  mass.  The  lower  bound  conditions 
given  below  in  this  section  are  not  violated.  Stars  flagged  in 
Hipparcos  with  “C”  (resolved  visual  binaries)  or  “O”  (pre¬ 
viously  known  astrometric  binaries  with  orbital  fits)  are  not  in¬ 
cluded  in  our  catalog,  but  a  small  number  of  stochastic  solutions 
(flag  “X”)  may  appear  there.  The  nearest  astrometric  binaries  of 
long  periods  are  of  special  interest,  since  they  may  have  substellar 
or  planetary  companions. 

Equation  (A13),  shown  in  Appendix  A,  can  be  rewritten  in 
terms  of  the  orbital  period  P, 


A  p  < 


2ttHRoM2 

Ml {V/3  ’ 


(1) 


where  n  is  the  parallax,  Mi  is  the  secondary  mass,  and  Mtot  is 
the  total  mass.  From  this,  given  a  period  estimate,  a  lower  es¬ 
timate  for  the  secondary  mass  can  be  derived.  The  upper  bound 
for  Ap  is  proportional  to  q2  =  M2I Mtot  and  has  a  weak  depen¬ 
dence  on  Mtot,  in  that  Apmax  oc  q2M^ .  The  inequality  be¬ 
comes  an  exact  equation  for  face-on  orbits.  The  projected  Ap  of 
an  inclined  orbit  is  always  smaller  than  the  true  center-of-mass 
deviation  of  proper  motion;  hence,  generally,  heavier  compo¬ 
nents  are  required  to  generate  the  same  astrometric  effect. 

Another  kind  of  uncertainty  is  related  to  the  fact  that  the  above 
formula  is  valid  for  the  difference  between  instantaneous  proper 
motion  and  the  true  center-of-mass  proper  motion,  whereas  the 
proper  motions  in  the  catalogs  are  some  average  of  proper  mo¬ 
tions,  over  about  3.2  yr  for  Hipparcos  and  several  decades  to  a 
century  for  Tycho-2.  If  the  period  of  a  system  is  longer  than  a  few 
decades,  the  Tycho-2  proper  motion  may  be  very  different  from 
the  true  proper  motion.  On  the  other  hand,  if  the  period  is  shorter 
than  several  years,  the  Hipparcos  proper  motion  is  an  average  of 
instantaneous  proper  motions  at  individual  times  of  observation. 
Such  an  average  will  be  between  the  maximum  and  the  minimum 
instantaneous  A/u-values,  as  described  by  the  Rq  (orbital  phase) 
factor  (Fig.  1).  The  averaging  effect  of  the  finite  duration  of 
Hipparcos  measurements  can  be  imagined  as  a  sliding  mean  of 
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Fig.  1. — Factors  R0,  Ri9  R2,  and  R2/R1  as  functions  of  orbital  phase  (t  —  T)/P  for  eccentricities  e  —  0  {dotted  lines),  e  =  0.2  {dashed  lines),  e  =  0.4  {dash-dotted 
lines),  and  e  —  0.6  {solid  lines). 


the  curves  in  Figure  1.  For  eccentric  orbits,  instantaneous  proper 
motion  deviations  are  greater  at  times  of  periastron  passages; 
therefore,  astrometric  binaries  of  high  eccentricity  can  be  better 
detected  with  this  method.  The  averaging  effect  will  diminish  the 
advantage  of  eccentric  binaries,  as  the  Rq  curves  will  become 
smoother  and  closer  to  unity.  Thus,  no  significant  overrepresen¬ 
tation  of  eccentric  A/i  binaries  is  expected  in  the  sample. 

For  very  long  period  systems,  the  Tycho-2  and  Hipparcos 
proper  motions  will  yield  a  A p  approaching  0  as  P  increases 
beyond  several  centuries.  This  case  was  analyzed  in  Kaplan  & 
Makarov  (2003),  and  it  was  shown  that  the  change  in  measured 
proper  motion  is  linear  with  the  observational  time  span.  The 
A/i  method  will  clearly  miss  those  systems,  as  the  apparent 
acceleration  will  be  too  small  to  produce  A/i  larger  than  2- 
3  mas  yr  1 .  Perhaps  the  next-generation  astrometric  satellites 
(SIM  and  Gaia)  will  detect  many  of  those  systems  at  the  mi- 
croarcsecond  level  of  precision  as  slowly  accelerating  stars. 

Another  limitation  of  the  astrometric  method  is  the  implicit 
assumption  that  the  secondary  component  is  much  dimmer  than 
the  primary  and  therefore  that  the  observable  photocenter  co¬ 
incides  with  the  primary,  as  discussed  in  the  appendices.  In 
some  cases,  when  the  brightness  ratio  can  be  determined  pho¬ 
tometrically  or  spectroscopically,  appropriate  corrections  can 
be  easily  introduced  into  the  formulae  relating  the  astrometric 
observables  and  the  orbital  elements.  In  the  marginal  case  in 
which  both  components  have  the  same  brightness,  the  photo¬ 
center  will  show  no  orbital  motion.  Most  companions,  however, 
have  significantly  smaller  masses  than  the  primaries,  and  since 
the  intensity  oflight  from  a  solar-type  star  with  Mg  [0.5,  2.0]  is 
approximately  proportional  to  M5-58  (Flenry  &  MacCarthy 
1993),  the  vast  majority  of  main-sequence  pairs  have  small 
photocentric  effects.  The  inequalities  given  in  this  section  and 


throughout  the  paper  are  valid  for  binaries  with  considerable 
photocentric  displacements,  because  they  specify  the  lower 
bounds  of  the  apparent  orbit  size  and  motion  and  the  secondary 
mass  ratio. 

Equation  (1)  shows  that  nearby  (II  >  40  mas)  solar-type  stars 
with  low-mass  companions  (</2  =  M^IM^  as  low  as  0.05-0.1) 
can  be  detected  as  Ap  binaries.  This  allows  the  possibility  of 
discovering  new  brown  dwarfs  around  nearby  stars  with  orbital 
periods  larger  than  4  yr,  still  a  poorly  known  domain  of  sub- 
stellar  companions  (Liu  et  al.  2002).  It  may  ultimately  reveal 
whether  the  “brown  dwarf  desert,”  the  apparent  paucity  of 
brown  dwarfs  with  short-period  orbits  around  solar-type  stars, 
extends  into  the  domain  of  longer  periods.  Companions  similar 
to  ice  giant  planets  can  be  detected  around  the  nearest  and  the 
smallest  stars,  such  as  M  dwarfs. 

Equation  (1)  is  readily  transformed  into 


<?2  > 


0o 

Ro 


(2) 


where  Qo  =  A/i/(27rII).  Here  the  factor  Q0  contains  the  two 
observables,  while  the  remainder  of  the  variables  are  unknown 
but  often  subject  to  reasonable  estimates  or  can  be  bounded. 
Decimal  logarithms  of  Qq  are  listed  in  the  catalog  of  A/i  bi¬ 
naries  (Table  1).  Whenever  the  parallax  was  less  than  2  mas,  we 
assumed  n  =  2  mas  and  gave  an  uncertainty  flag  “U”  in  the 
end  of  the  record.  Those  distant  stars,  interestingly,  have  large 
lower  bounds  of  q2,  implying  either  that  the  dark  companion 
is  actually  more  massive  than  the  primary  or  that  the  stars  are 
quite  eccentric,  and  that  the  observation  took  place  near  the 
periastron  time.  In  each  particular  case  of  interest,  a  more  careful 
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study  should  be  done,  since  other  occasional  perturbations  can¬ 
not  be  precluded.  For  example,  a  small  fraction  of  distant  stars 
may  be  unresolved  visual  binaries  with  a  slowly  variable  com¬ 
panion.  This  may  cause  the  photocenter  to  shift  slowly  along  the 
line  connecting  the  components  (the  effect  known  as  a  variability- 
imposed  mover),  producing  apparent  proper  motion  or  even  ac¬ 
celeration.  The  majority  of  objects  in  the  catalog  are  genuine 
binaries,  in  some  cases  with  low-mass,  dim  companions. 

Consider,  as  an  example,  the  star  AB  Dor  (HIP  25647).  It 
is  one  of  the  very  young  and  active  stars  in  the  immediate  so¬ 
lar  neighborhood,  whose  origin  is  unknown  (  Makarov  2003). 
Guirado  et  al.  (1997)  combined  VLBI  and  Hipparcos  data  to 
derive  a  putative  orbit  for  this  astrometric  binary.  According  to 
their  estimation,  the  orbital  period  is  between  6.5  and  27.5  yr. 
From  the  catalog  of  Ag  binaries,  we  have  log  Qq  =  —1.3.  As¬ 
suming  a  total  mass  ofO.86M0,  we  obtain  q2  >  0.10/RoatP  = 
6.5,  and  q2  >  0.16/7?o  at  P  =  27.5  yr.  The  analysis  of  Guirado 
et  al.  implies  that  the  Hipparcos  observations  were  taken  around 
the  periastron  time,  and  the  eccentricity  is  between  0.28  and 
0.78.  Hence,  the  orbital  phase  factor  R0  can  be  between  1.3  and2.8. 
Taking  this  factor  into  account  reverses  the  boundaries  of 
the  possible  interval,  so  that  <72  >  0.07  at  P  —  6.5,  and  q2  > 
0.06  at  P  =  27.5  yr.  We  finally  conclude  that  the  secondary 
mass  in  this  system  can  be  as  small  as  0.05-0.06  M0,  in  good 
agreement  with  the  estimates  in  Guirado  et  al.  (1997)  of  0.08- 
0.11.  The  actual  mass  may  indeed  be  somewhat  larger  than  our 
lower  limits,  since  the  inclination  of  the  orbit  i  is  probably 
about  60°. 


3.  fi  BINARIES  AND  fi  BINARIES 

All  currently  known  g  and  fl  binaries  are  collected  in  a  cat¬ 
alog  of  2622  stars  (Table  2).  The  data  are  mostly  copied  from 
the  Double  and  Multiple  Star  Annex,  Part  G,  of  the  Hipparcos 
catalog.  We  added  only  cross-reference  flags  “M”  for  stars  that 
are  also  Ag  binaries,  estimates  related  to  the  mass  ratio  lower 
bounds,  and  uncertainty  flags  “U”  for  the  latter  whenever  the 
Hipparcos  parallax  was  less  than  2  mas.  It  is  noted  again  that 
the  orbital  estimates  (Appendix  B)  suffer  three  types  of  un¬ 
certainties:  (1)  the  projection  effect  for  inclined  orbits,  (2)  the 
averaging  effect  due  to  the  finite  duration  of  the  observation, 
and  (3)  the  orbital  phase  effect  for  eccentric  orbits.  The  former 
effect  always  makes  the  observed  accelerations  smaller  than 
the  true  orbital  accelerations,  and  thus  our  upper  estimates  of  a 
larger.  The  latter  two  effects  may  act  either  way.  In  order  to 
better  illustrate  these  effects,  the  orbital  phase  factors  Rq,  R\, 
and  R2  are  depicted  in  Figure  1  as  functions  of  phase  for  e  =  0. 
0.2,  0.4,  and  0.6.  The  heights  of  the  peaks  near  the  periastron 
increase  for  accelerating  binaries  and  stars  with  significant 
second  derivatives;  hence,  the  effect  of  eccentricity  becomes 
stronger.  Although  the  duration  of  enhanced  nonlinear  motion 
around  the  periastron  time  is  only  about  1/10  of  the  orbital 
period,  eccentric  binaries  that  happen  to  be  observed  close  to 
their  periastron  have  much  better  chances  of  being  detected. 
Indeed,  the  magnitudes  of  g  and  fi  are  proportional  to  the  par¬ 
allax  n  (eqs.  [B12]  and  [B 1 7]),  while  the  number  of  detectable 
eccentric  binaries  grows  as  n~3  in  the  near-solar  neighborhood. 
Thus,  binaries  of  extreme  eccentricity  could  be  detected  even  at 
relatively  large  distances. 

Equation  (B14)  from  Appendix  B  can  be  rewritten  in  terms 
of  the  orbital  period  P\ 


.  (2ii?i)2nM2 


(3) 


From  this,  given  a  period  estimate,  a  lower  estimate  for  the  secon¬ 
dary  mass  can  be  derived.  Just  as  for  Ag,  the  upper  bound  for  fi  is 
proportional  to  q2  =  M2/ Mtot  and  has  a  weak  dependence  on  Mtot , 
in  that  /imax  oc  q2M^ .  The  inequality  becomes  an  exact  equation 
for  face-on  orbits.  Equation  (3)  is  readily  transformed  into 


qi  > 


Qi  P4/3 
Rj  M1/3 

1  mtot 


(4) 


where  Q\  —  /i/(47r2n).  Decimal  logarithms  of  (Q\  are  listed  in 
the  catalog  of  fi  binaries  ( Table  2).  Whenever  the  parallax  was 
below  2  mas,  we  assumed  n  =  2  mas  and  gave  an  uncertainty 
flag  “U”  in  the  end  of  the  record. 

Figure  1  ( bottom  right)  shows  why  second  derivatives  of 
proper  motions  can  be  larger  in  magnitude  than  accelerations, 
even  in  long-period  binaries.  This  is  caused  by  the  large  cur¬ 
vature  of  highly  eccentric  orbits  near  the  periastron.  Since  the 
ratio  R2/R\  is  significantly  larger  than  unity  only  within  ±0.  IP, 
the  probability  of  observing  such  accelerations  is  fairly  small 
for  a  given  eccentric  binary.  Presumably,  very  eccentric  binaries 
are  rare,  too.  The  fraction  of  stars  with  jil g  >  1  should  be  quite 
small  in  an  unbiased  sample.  However,  most  of  the  stars  with 
significant  accelerations  and  second  derivatives  have  ft/ g  >  1, 
which  is  probably  a  statistical  bias  due  to  the  much  higher  de¬ 
tection  threshold  for  jl.  Stars  with  significant  fi  can  hardly  have 
periods  P  »  1  yr,  as  discussed  in  §  4. 

The  statistical  significances  of  accelerations  and  second  deriv¬ 
atives  are  given  as  statistics  in  columns  (7)  and  (10)  of  Table  2, 
respectively.  These  statistics  are  equivalents  to  \2  statistics  with 
2  degrees  of  freedom  for  the  case  of  a  multivariate  normal  dis¬ 
tribution.  Estimated  accelerations  were  considered  statistically 
significant  when  that  quantity  was  larger  than  3.44,  and  similarly 
for  second  derivatives  of  proper  motion.  This  criterion  corre¬ 
sponds  to  the  3  a  level  of  significance  of  a  one-dimensional 
Gaussian  distribution,  or  a  0.27%  probability  of  a  null  hypothe¬ 
sis  (that  the  true  parameter  value  is  0). 


4.  OVERVIEW  OF  ASTROMETRIC  BINARIES 
From  relations  (A13)  and  (B12),  one  obtains 

—  (1  +  ecosP)^2(l  —  ecosP)3//2.  (5) 

p  2tt 

This  equation  is  strict  for  face-on  orbits  and  a  good  approxi¬ 
mation  for  moderately  inclined  orbits.  For  circular  orbits,  the 
ratio  of  instantaneous  proper  motion  difference  and  acceleration 
is  roughly  P/2tt.  This  approximate  equality  is  expected  to  hold 
for  the  majority  of  astrometric  binaries  of  moderate  eccentric¬ 
ities,  since  the  probability  of  catching  a  star  near  its  periastron  is 
small.  However,  very  eccentric  binaries  may  tip  the  balance 
toward  stars  with  small  Ag/ fi.  With  that  in  mind,  we  discrim¬ 
inate  in  the  set  of  astrometric  binaries  in  the  two  catalogs  a  set  of 
binaries  with  significant  Ag  but  nonmeasurable  g,  and  a  set  of 
binaries  with  significant  g  but  insignificant  A g.  It  is  expected 
that  most  of  the  former  should  have  periods  longer  than  about 
6  yr,  while  the  latter  group  should  include  mostly  stars  with  pe¬ 
riods  shorter  than  6  yr. 

Figure  2  shows  distributions  of  B  —  V  colors  of  the  short- 
period  and  long-period  binaries,  selected  as  described  above,  and 
the  general  distribution  of  Hipparcos  stars  ( solid  line )  normal¬ 
ized  to  unit  area.  There  is  no  appreciable  difference  between  the 
distributions  of  short-period  binaries  and  the  general  sample. 
Hence,  short-period  binaries  have  the  same  rate  of  appearance  in 


TABLE  2 


A  Catalog  of  fi  Binaries  and  p  Binaries 


HIP  No. 

(1) 

TYC1 

(2) 

TYC2 

(3) 

TYC3 

(4) 

fta* 

(5) 

Vs 

(6) 

Sig.  Level 
(7) 

Ma* 

(8) 

Vs 

(9) 

Sig.  Level 
(10) 

n 

(ID 

A/i 

Binary? 

(12) 

log  Q\ 
(13) 

Uncertain? 

(14) 

62 . 

8022 

411 

1 

-6.41  ±2.69 

-5.16  ±  1.71 

3.53 

5.3  ±  1.1 

-1.3 

68 . 

1178 

1142 

1 

3.04  ±  2.50 

-7.04  ±  1.66 

5.33 

31.8  ±  1.2 

M 

-2.5 

104 . 

6989 

95 

1 

9.75  ±2.79 

5.97  ±  1.38 

5.63 

8.2  ±  1.2 

-1.3 

305 . 

6418 

1218 

1 

-4.26  ±  2.97 

9.95  ±  1.40 

7.11 

20.4  ±  1.0 

M 

-2.2 

340 . 

7529 

819 

1 

3.16  ±  3.19 

-0.99  ±  2.99 

1.15 

-29.08  ±  8.87 

-17.63  ±  6.97 

3.84 

14.9  ±  1.0 

-2.4 

493 . 

1181 

1071 

1 

-0.60  ±  1.63 

3.96  ±  1.34 

3.58 

26.2  ±  0.8 

-2.5 

611 . 

7526 

641 

1 

-17.65  ±  3.33 

-2.38  ±  3.24 

5.42 

19.19  ±  9.25 

-25.18  ±  7.91 

4.15 

7.2  ±  1.1 

-1.2 

646 . 

5841 

423 

1 

-7.51  ±2.53 

-2.83  ±  1.02 

4.02 

3.6  ±  0.9 

-1.1 

648 . 

9350 

1236 

1 

6.77  ±  1.80 

3.36  ±  1.52 

4.19 

-17.15  ±  4.31 

9.84  ±  4.40 

4.96 

8.9  ±  0.6 

-1.5 

695 . 

4018 

3889 

1 

0.18  ±  1.13 

6.17  ±  1.66 

3.73 

2.2  ±  0.6 

-1.1 

741 . 

6995 

264 

1 

11.84  ±  3.24 

5.29  ±2.38 

4.07 

0.1  ±  1.4 

-0.4 

U 

Notes. — Several  entries  are  shown  as  an  example;  the  complete  version  of  the  catalog  is  available  through  the  CDS  Web  site.  Col.  (1):  Hipparcos  numbers.  Cols.  (2) — (4):  Tycho  identification  numbers  TYC1,  TYC2, 
and  TYC3.  Cols.  (5)  and  (6):  Hipparcos  acceleration  components  paif  and  ps  and  their  errors  in  mas  yr-2.  Col.  (7):  Significance  levels  of  acceleration  terms.  Cols.  (8)  and  (9):  Second  derivatives  of  proper  motion  paie  and 
p$  and  their  errors  in  mas  yr-3.  Col.  (10):  Significance  levels  of  second  derivatives.  Col.  (1 1):  Hipparcos  parallaxes  and  their  errors  in  mas.  Col.  (12):  Flag  “M”  is  set  for  stars  that  are  also  A/i  binaries.  Col.  (13):  log  Q\ 
factors  for  the  lower  limits  of  q 2.  Col.  (14):  Flag  “U”  for  uncertain  log  Q\  estimates  for  stars  with  small  or  negative  parallaxes. 
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B-V  mag 

Fig.  2. — Distributions  of  B  —  V  colors  for  all  1 1 8,2 1 8  Hipparcos  stars  (solid 
line),  astrometric  binaries  with  significant  Ap  differences  but  insignificant  ac¬ 
celerations  (dot-dashed  line;  1161  stars),  and  with  significant  accelerations  but 
insignificant  Ap  (dashed  line;  1854  stars). 

all  spectral  types.  This  is  obviously  not  the  case  with  long-period 
astrometric  binaries,  which  tend  to  appear  mostly  as  solar-type 
spectral  types  (late  F  and  G)  and  rarely  as  B-  and  A-type  stars. 
It  may  be  tempting  to  interpret  this  result  in  terms  of  different 
binary  formation  mechanisms  (e.g.,  fragmentation  vs.  capture) 
producing  distinct  period  distributions  for  more  massive  stars. 
We  give  a  cautionary  note  that  early-type  stars  are  also  younger, 
and  some  evolutionary  effects  cannot  be  precluded. 


of  stars  in  which  a  reasonable  degree  of  completeness  can  be 
anticipated.  Some  aspects  of  binary  formation  at  low  q2  can  be 
addressed,  such  as  distributions  of  orbital  elements  and  their 
dependencies  on  age  and  metallicity. 

The  fraction  of  astrometric  binaries  will  dramatically  increase 
when  the  next-generation  astrometry  projects  provide  us  with 
much  more  accurate  data  ( Gaia  and  SIM).  In  fact,  astrometric 
binaries  may  prove  to  be  a  considerable  difficulty  in  the  com¬ 
plicated  data  reduction  systems  for  these  missions.  If,  for  ex¬ 
ample,  brown  dwarf  companions  are  widespread  in  long-period 
binaries,  a  large  number  of  reference  stars  may  be  completely 
unsuitable  for  processing  with  the  standard  astrometric  model  of 
linear  motion  and  parallax. 


The  research  described  in  this  paper  was  carried  out  in  part  at 
the  Jet  Propulsion  Laboratory,  California  Institute  of  Technol¬ 
ogy,  under  a  contract  with  the  National  Aeronautics  and  Space 
Administration. 

APPENDIX  A 

ORBIT  CONSTRAINTS  FROM  A p 

The  apparent  motion  of  a  binary  in  the  plane  of  celestial 
projection  is  described  by  (Heintz  1978) 

x  =  A(cosE  —  e)  +  FV 1  —  e2  sin  is,  (Al) 

y  —  B(cosE  —  e)  +  GV 1  —  e2  sin  is,  (A2) 


5.  SUMMARY 

The  two  catalogs  of  astrometric  binaries  presented  here,  based 
on  the  results  of  the  Hipparcos  mission  and  older  ground-based 
astrometric  catalogs,  can  be  useful  in  the  ongoing  quest  for 
low-mass  binaries  and  brown  dwarfs  in  the  solar  neighbor¬ 
hood.  The  method  of  astrometric  motion  analysis  is  sensitive  in 
the  most  difficult  area  of  orbital  parameters  for  both  spectro¬ 
scopic  and  imaging  investigation,  i.e.,  orbital  periods  between 
3  and  100  yr.  Although  the  proper  motion  data  are  not  sufficient 
to  estimate  the  orbital  elements,  approximate  constraints  are 
useful  for  statistical  studies  and  for  selecting  promising  candi¬ 
dates  among  nearby  stars.  In  most  cases,  additional  information 
via,  e.g.,  radial  velocity  monitoring  or  new  astrometric  mea¬ 
surements  will  be  required  to  prove  the  existence  of  a  low-mass 
companion  and  to  estimate  its  mass.  The  nearest  binaries  are 
especially  interesting,  since  the  accuracy  of  the  Hipparcos  and 
Tycho-2  data  is  sufficient  to  detect  giant  planets.  However,  care 
should  be  exercised  with  nearby  stars  that  have  large  radial 
velocities,  since  their  large  parallaxes  and  considerable  secular 
accelerations  were  not  taken  into  account  in  deriving  the  Tycho-2 
proper  motions.  For  example,  existence  of  planetary  companions 
to  Barnard’s  star  and  Kapteyn’s  star  should  not  be  taken  for 
granted  just  because  these  stars  are  included  in  the  catalog  of 
astrometric  binaries.  Barnard’s  star,  in  particular,  has  a  large  sec¬ 
ular  acceleration  of  1.24  mas  yr~2  in  declination  that  amounts 
to  an  astrometric  offset  of  6".2  over  100  yr.  This  offset  may 
account  for  the  «50  mas  yr~*  difference  between  the  Tycho-2 
and  Hipparcos  proper  motions.  Such  objects  require  more  so¬ 
phisticated  astrometric  analysis. 

Because  of  the  limited  accuracy  of  the  available  astrometric 
data,  the  number  of  astrometric  binaries  is  currently  not  large. 
However,  important  conclusions  can  be  drawn  from  the  avail¬ 
able  sample,  especially  on  a  well-defined  volume-limited  sample 


where  x  and  y  are  the  coordinates,  e  is  the  eccentricity  of  the 
orbit,  and  E  is  the  eccentric  anomaly  related  to  the  mean  anom¬ 
aly  M  by  Kepler’s  equation 

t  -  T 

M  =  It: — — —  =  E  —  esinE.  (A3) 


As  usual,  T  is  the  periastron  epoch  and  P  is  the  orbital  period. 
The  Thiele-Innes  constants  are  related  to  the  remaining  orbital 
elements  by 

A  =  aa( cos  lu  cos  Q  —  sin  lu  sin  f l  cos  i ),  (A4) 

B  =  aa( cos  lu  sin  f2  +  sin  lu  cos  H  cos  i),  (A5) 

F  =  aa(—  sin  w  cos  fl  —  cos  w  sin  H  cos  i),  (A6) 

G  =  aa(—  sin  cu  sin  fi  +  cos  lu  cos  fl  cos  i),  (A7) 


where  aa  is  the  astrometric  semimajor  axis  (whether  for  the 
primary,  secondary,  or  the  photocenter),  lu  is  the  longitude  of 
the  periastron  in  the  plane  of  orbit,  and  O  is  the  position  angle 
of  the  node  in  the  plane  of  projection. 

The  first  derivative  of  E  is 


E  = 


2n  1 
P  1  —  ecos  E 


(A8) 


Differentiating  equation  (A2)  with  respect  to  time  obtains  (see 
also  Brouwer  &  Clemence  1961) 


x  =  (—A  sin  E  +  FV  1  —  e2  cos  E)  E ,  (A9) 

y  =  (-B  sin E+  GVl  -  e2  cos  E)E.  (A10) 
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The  square  of  the  apparent  instantaneous  orbital  velocity  is 

,  2 


x2  +  y2  =  a2a 


^cos  u)  sin  E  +  sin  ujV  1  —  e2  cos  E^j 
+  ^sin  to  sin  E  —  cos  uV  1  —  e2  cos  E^j  cos2/ 


E2. 


(All) 

Note  that  the  magnitude  of  the  instantaneous  velocity  is  inde¬ 
pendent  of  12.  Since  cos2/  <  1,  and  for  the  primary  component 


aa  =  a\U.  = 


aHM2 

V/fnt 


(A12) 


we  obtain  after  replacing  cos 2i  with  unity. 


and  a  is  the  semimajor  axis  of  the  relative  orbit.  This  can  be  re¬ 
arranged  to 


2  f/'gr  -Hot  /  2.U  \ 

a  \r  ) 


(A18) 


where  r  =  a(  1  —  e2)/(l  +  e  cos  9)  and  9  is  the  angle  at  the  pri¬ 
mary  between  the  direction  to  periastron  and  the  instantaneous 
position  of  the  secondary;  i.e.,  9  is  the  true  anomaly.  In  terms  of 
eccentric  anomaly,  r  =  a(l  —  e  cos  E ).  It  is  easy  to  show  that 
the  factor  in  parentheses  above  equals  R jj. 

For  face-on  orbits,  the  projection  factor  is  1  everywhere,  and 
A fi  =  vHM2/Mtot.  Therefore,  we  have 


/GgrMpA  1/2HR0M2 

\  a  J  dTt0t 


(A19) 


27rII  R0M2 

Ap,  <  — 7^=- , 

v  oMtot 


(A13)  and  if  GgI  =  4i r2,  this  becomes 


where  II  is  the  parallax  in  the  same  units  as  the  semimajor  axis 
of  the  primary  aa,  a  1  is  the  semimajor  axis  of  the  primary  in  AU, 
M2  is  the  mass  of  the  secondary,  Mtot  is  the  total  mass,  A/x  is 
the  astrometrically  measured  instantaneous  orbital  velocity  of 
the  primary  in  the  plane  of  projection,  with  respect  to  the  center 
of  mass,  and  the  time-dependent  (orbital  phase)  term  is 


Ap 


2ttHRqM2 
v2  ci  Mtot 


(A20) 


the  same  as  equation  (A  13)  applied  to  face-on  orbits. 
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1  +  ecos  E\  1  /  2 
1  —  e  cos  E ) 


(A14) 


This  can  be  transformed  into  a  constraint  for  the  full  semimajor 
axis  of  the  orbit: 


z  (27r)2n2R2M2 

0  ~  MotA/i2 


(A15) 


If  the  share  of  intensity  of  the  secondary  star  in  the  total  in¬ 
tensity  of  the  system,  k2  =  I2  /  (I\  + 12),  is  not  negligible  and  the 
separation  between  the  components  is  much  smaller  than  the 
resolution  limit  of  the  instalment,  M2  in  equations  (A13)  and 
(A15)  should  be  replaced  with  {k\M2  —  k2M\).  In  this  case,  the 
astrometric  excursion  of  the  photocenter  is  always  smaller  than 
the  actual  excursion  of  the  primary,  and  subsequently,  the  ob¬ 
served  A/x  is  smaller  too. 

For  Hipparcos  stars  that  have  significant  differences  in  proper 
motion  between  the  long-term  Tycho-2  and  Sixth  Fundamental 
Catalogue  but  do  not  have  measurable  accelerations,  a  lower 
bound  of  a  can  be  derived  from  the  somewhat  arbitrary  condi¬ 
tion  that  the  orbital  period  should  be  greater  than  or  approxi¬ 
mately  equal  to  5  times  the  mission  duration  (typically,  3.2  yr). 
This  provides  a  rough  condition: 


ORBIT  CONSTRAINTS  FROM  /2 
Differentiating  equation  (A10)  with  respect  to  time  obtains 

x  =  —  (^A  cos  E  +  FV 1  —  e2  sin  eJ  E2 

+  (—A  sin  E  +  FV l-e2  cos E'jE,  (Bl) 

y  =  —  (^B  cos E  +  GV 1  —  e2  sin E^j E 2 

+  B  sin  E  +  GV l-e2  cos  Ej  E.  (B2) 


Substituting 


•  ?  1  —  e  cos  E  ■■ 

E2  = - - - E. 


e  sin  if 

squaring,  and  adding  up  equations  (Bl)  and  (B2)  obtains 


(B3) 


T7  2 


fi2  =x 2  +y2  =  —  (piq2  -  2p2q\q2  +  p2q f), 


(B4) 


where 


a>63M[J2.  (A16) 

An  alternative  way  to  derive  the  equality  corresponding  to 
equation  (A13)  is  to  use  the  vis  viva  integral,  representing  con¬ 
servation  of  energy, 

v2  =  2GgrMm^--^-]j,  (A17) 

in  which  r  is  the  instantaneous  distance  between  the  two  com¬ 
ponents,  v  is  the  relative  speed,  Ggr  is  the  gravitational  constant. 


pi  =  A2  +  B2  =  a2  (cos  2lo  +  sin2u)  cos2/) ,  (B5) 

p2  =  F2  +  G2  =  a2  (sin2w  +  cos2w  cos2/) ,  (B6) 

p2  =  AF  +  BG  =  a2  (—cos  tu  sin  lu  +  cos  u>  sin  oj  cos2/) , 


9\  = 
<72 


1  —  cos  E 
sin£  ’ 

=  Vl  —  e2. 


(B7) 

(B8) 

(B9) 
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Substituting  the  expressions  for  pu  p2,  and  p2  leads  to  a  qua¬ 
dratic  form  similar  to  equation  (All): 


x2  +  v2  =  a^rE1  (cos  ujq\  +  sin  u>q2)2 
eL  L 

+  (sincucjq  —  cos  u>q2)2  cos2/ 


(BIO) 


Since  cos2/  <  1,  using  equation  (A12)  and 


For  Hipparcos  stars  that  have  significant  proper  motion  ac¬ 
celerations  but  do  not  have  orbital  solutions,  a  lower  bound  of  a 
can  be  derived  from  the  somewhat  arbitrary  condition  that  the 
orbital  period  should  be  greater  than  or  approximately  equal  to 
2  times  the  mission  duration  (typically,  3.2  yr).  This  provides  a 
rough  condition: 


a>3AM^3. 


(B15) 


e  sin  E 


P  )  (1  -  e  cos  Ef  ' 


(Bll) 


we  derive  for  the  instantaneous  acceleration  of  proper  motion 
in  the  plane  of  projection 


R 


i  =  \Jx2  +  y2  < 


/  2nR\ 

\  a 


UM2, 


(B12) 


where  II  is  the  parallax  in  the  same  units  as  the  astrometric 
excursion  aa,  M2  is  the  mass  of  the  secondary  component,  and 
the  time-dependent  term 


Ri  = 


1 

1  —  ecos  E 


(B13) 


Alternatively,  the  equality  relation  (B 12),  corresponding  to  a 
face-on  orbit  case,  can  be  derived  directly  from  the  formula  for 
orbital  acceleration, 


Ggr  A/tot 

r2  ’ 


(B16) 


along  the  lines  described  at  the  end  of  Appendix  A,  using  the 
relation  fi  =  i)HM2/Mtot  for  instantaneous  acceleration  of  the 
primary  with  respect  to  the  center  of  mass. 

We  omit  the  somewhat  cumbersome  derivation  for  fl  and 
give  the  final  formula: 


l'j  = 


fdlx\2  fdW  ^  dn)2 UM2^M^tR2 
\dt2)  \dt2J  ~  aP 2 


(B 17) 


This  can  be  transformed  into  a  constraint  for  the  full  semimajor 
axis  of  the  orbit: 


a  <  2ttR\ 


(B14) 


where  the  time-dependent  term 

\J  1  +  e2(3  —  4  cos 2E) 
(1-ecosis)4 


(B18) 


If  the  share  of  intensity  of  the  secondary  star  in  the  total  in¬ 
tensity  of  the  system,  k2  =  I2I  (I\  + 12),  is  not  negligible  and  the 
separation  between  the  components  is  much  smaller  than  the 
resolution  limit  of  the  instrument,  M2  in  equations  ( B 12)  and 
(B14)  should  be  replaced  with  ( k\M2  —  k2M\).  In  this  case,  the 
astrometric  excursion  of  the  photocenter  is  always  smaller  than 
the  actual  excursion  of  the  primary,  and  subsequently,  the  ob¬ 
served  [i  is  smaller  too. 


This  can  be  transformed  into  a  constraint  for  the  full  semimajor 
axis  of  the  orbit: 

^  (27r)6/7(R2nM2)2/7Mt10{7 
«< - ^2/7 - •  (B19) 

Inequalities  (B 1 7)  and  (B19)  become  exact  equations  for  face- 
on  orbits  (/  =  0°). 
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